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The 2d and 3d like Gletzer, Okhitani and Yamada (GOY) shell models are examined. The 2d 
like model shows a transition from statistical quasi-equilibrium to cascade of enstrophy as a function 
of the spectral ratio of energy to enstrophy. The transition is related to the ratio of time scales, 
corresponding to eddy turnover times, between shells. The anomalous scaling, giving rise to non- 
linear scaling functions, is also connected to the ratio of eddy turnover times. This is illustrated 
in a simple stochastic model, where the structure function, £(q), becomes independent of q. In the 
3d like model the multiscaling is also influenced by the existence of a second non-positive definite 
ON ' inviscid invariant, the helicity. 

The problem of understanding the scaling properties of velocity correlations in isotropic and homogeneous turbulence 
is still largely unsolved. The only exact derived property is the Kolmogorov scaling law for the third order moments, 
O . Cs(r) — —4er/5 + 6vdC2(r)/dr, where C q (r) = (Av(r) q ) = (\v(x + r) — v(r)\ q ),r is the distance in the fluid and (.) is 
the ensamble average. This relation is not closed, since it contains both second and third order correlation functions. 
In the inertial range defined by, L r ^ A, where L is the outer scale and A the inner or Kolmogorov dissipation 
£NJ ' scale, the second term on the right hand side is negligible and the scaling relation, C%{r) ~ r^ 3 ', with £(3) = 1, holds. 
The classical Kolmogorov theory (K41) based on dimensional arguments states that C q (r) ~ r^^ q \ with £(q) = q/3. 
It was noted by Landau shortly after the theory was presented that the energy dissipation could vary so much as to 
alter the scaling laws. This was incorporated in a refined version of Kolmogorovs theory (K62), in which a dependence 
' on L/r, the ratio of the outer scale to the scale in the inertial range was incorporated. 

It has been seen in numerous experiments that £(q) is a weakly non-linear function of q, different from the Kol- 
mogorov predictions for q different from 3. It follows from the Holder inequality that is a convex of q, thus 
> q/3 for q < 3 and ((q) < q/3 for q > 3. The deviation of the scaling from the K41 prediction is referred to as 
intermittency corrections. It is widely attributed to the fact that the energy dissipation in fully developed turbulence 
is indeed highly inhomogeneous, basically taking place on lower dimensional subsets, filaments, of the flow field. 

Numerical studies of the GOY shell model jl| of turbulence have been popular resently, mainly because this model, 
in which scaling relations completely equivalent to those of turbulence, also shows intermittency corrections to the 
K41 - equivalent - predictions. The GOY model structually resembles the spectral form of the Navier-Stokes equation, 
but there is no spatial fields associated with the wave components in the GOY model. The intermittency corrections 
in the scaling properties are thus associated with temporal intermittency which can only be weakly linked to spatial 
, intermittency through the Taylor hypothesis. The connection between the spatial intermittency of turbulence and 
temporal intermittency of the GOY model is not clear, but the hope is that understanding the latter can shed light 
• i-H , upon the former. This note is about the latter. 

In the GOY model the spectral domain is represented as shells, each of which is defined by a wavenumber k n = A™, 
where A is a scaling parameter defining the shell spacing. There are 2N degrees of freedom, where N is the number 
of shells, namely the generalized complex shell velocities, u n for n — 1,N. The dynamical equation for the shell 
velocities is, 

u n = ik n (u n+ 2U n+ i - ^it n+ iii„-i + - ^ 2 ^ u n -iu n - 2 )* - vk 2 n u n + fS„ t n , (1) 

where the first term represents the non-linear wave interaction or advection, the second term is the dissipation, and 
the third term the forcing, where no is some small wavenumber. The boundary conditions are, u_i = uq = tijv+i = 
un+2 = 0. The model has two quadratic inviscid invariants, which are constants of motion in the case / = v = 0, 
E = Y^n \ u- n\ 2 , referred to as the energy, and H = X^i( e — l)~"|un| 2 - For e < 1 H is non-positive definite, referred to 
as the helicity, and the model is thought of as modeling 3d turbulence. For e > 1 H is positive definite, referred to 
as the enstrophy, and the model is meant to resemble 2d turbulence. Dimensional arguments similar to those of K41 
can be applied to the GOY model assuming dissipation of one of the conserved quantities. For the 2d type model, 
1 < e < 2, enstrophy will be dissipated, and an additional large scale drag term, —v'k~ 2 u ni is added to (1) in order to 
remove the energy. The reason for the enstrophy and not the energy being dissipated is the usual that the dissipation 
of energy is k~ a times the dissipation of enstrophy at shell n, thus negligible for a > and n — * oo. The Kolmogorov 
scaling for the 2d type model then becomes, (\u n \) ~ fc I 7*' Q+1 ' ) ^ 3 . Note that this is an unstable fixed point of (1) for 
/ = v = 0, describing a cascade. As pointed out by Aurell et al. B in the case of a = 2 this scaling coincides with 
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the scaling that would be obtained in a statistical equilibrium, in which the enstrophy is distributed evenly among 
the degrees of freedom of the system, (\u n \) ~ fc,7 Q ^ 2 . The a — 2, corresponding to e = 5/4, case is a borderline case 
between models showing cascade, for a < 2, and models showing statistical (quasi-)equilibrium, for a > 2. In order to 
show this a series of numerical model runs for different values of e have been performed, details of which are reported 
elsewhere ||. 
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FIG. 1. Results of runs of the 2d like GOY model for different values of e. Shown is the spectral slope 
7 = — log((it))/ log(fc) as function of a = — log(e — l)/log(A). The lines are the Kolmogorov scaling for enstro- 
phy cascade and equipartitioning of enstrophy, respectively. For the 2d like GOY models an additional dissipation 
of the form, —v'k„ 2 u n , is applied in order to remove energy at small wave number shells. The model is run with 
n = 50,fc = A" 4 , A = 2,/„ = 5xl(n 4 x(l+i)<5n,i5, (e = 11/10,^ = 5xl0~ 27 ,i/' = 100), (e = 117/100, v = 5xlO~ 2 V = 100), (e = 5/4, v = I 
The models show a cross-over at a = 2 between statistical equilibrium, a > 2, and enstrophy cascading, a < 2. 

Figure 1 shows the spectral slopes obtained. The two lines shows the slopes corresponding to cascade and statistical 
equilibrium respectively, the diamonds are the numerical results. The figure clearly shows a cross-over from one type 
of behavior to the other. The reason for this transition is related to the scaling of typical timescales, or eddy turnover 
times, for the different shells. The eddy turnover time at shell n is given as t„ = (kndunl))^ 1 as seen from (1) 
or dimensional arguments. The eddy turnover time then becomes, ^"cascade ^ 

respectively. Thus for a > 2 the eddy turnover times in both cases increase with wave number, and equilibration via 
inverse enstrophy transfer, from large wave numbered shells to small wave numbered shells takes place. For a < 2 the 
situation is reversed, and a statistical equilibrium can newer be achieved. Similar results have been found by Yamada 
and Ohkitani j| for a slightly different set of GOY models, with only one positive definite inviscid invariant, reffered 
here to as 3D type models. 

In the 3d type models the situation is totally different, in this case the second inviscid invariant, the helicity, is 
not cascaded even though the ratio of absolute value of the helicity to the energy at shell n grows exponentially 
with n as is the case for the enstrophy. The reason is that the helicity has alternating signs for even and odd 
numbered shells. Therefore there will be a net production of (positive sign) helicity at odd numbered shells due to 
the dissipation, and no net flow of helicity from small wave number shells to large wave number shells is necessary. It 
has been demonstrated numerically that helicity is indeed the quantity to be cascaded in the (pathological) case of 
hyperviscosity only active on the outermost shell |^] . In the usual 3d case energy will be cascaded, with the resulting 
spectral scaling, (|u| 3 ) ~ k~ x . 

Both the 2d and the 3d models shows intermittency corrections to the Kolmogorov scaling depending on e and A. 
In this note I will suggest two different mechanisms in play in the 3d case, where only the one is in action in the 2d 
case. 

In the 3d case it is observed that the structure function depends only on e and A in the combination, — log(l — 
e)/log(A) = a ||. With S((q) = C(<?) — <?/3, S((q) increases in absolute value when a increases. I suggest that this 
is due to differences in the ratio of helicity production and helicity elimination by dissipation at neighboring shells 
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in the beginning of the dissipative subrange where scaling still approximately holds or where extended self-similarity 
can be applied fjj. This ratio, r, can be estimated as r = (Aff n _i + AH n+1 )/2AH n ~ (A~ Q+2/3 + A Q ~ 2/3 )/2, 
which for (a, A) = (0.5, 2), (1, 2), (2, 2) gives 1.007,1.03,1.46 respectively. So there is the biggest "mis-match", or 
non-cancellation, in the case (2, 2) where the biggest non-linearity in the structure function is observed. That these 
two things should be related is consistent with the findings of a moderated GOY model (model 3) by Benzi et al. [p|, 
where 2 copies of the GOY model are coupled. In this model the helicity takes the form, H = J2 n ^n(l u n I 2 — \ u n I )> 
where u+ and u~ are the two complex variables in shell n. In this model the helicity production and elimination will 
- on average - exactly balance, thus no dependence of the structure function on a should be expected in agreement 
with the numerical findings. 

In the standard 3d case, a = 1, the structure function is still non- linear even though the above defined r is close 
to 1. This is suggested to be due to the difference in time scales between the large and the small scales. This effect 
is expected also to determine the anomalous scaling in the 2d case, where only the dissipation of enstrophy plays a 
role. In the 2d case numerical studies suggests that the absolute size of S((q) increases as a decreases from 2 to [[|. 
This suggests that intermittency is enhanced with growing difference between eddy turnover times from one shell to 
the next, which goes as t„ + i/t„ = A' 2- "^ 3 . 
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FIG. 2. The non-linear transfer of energy from one shell to the next in the 3d like GOY model with e = 1/2, all parameter 
values are standard, forcing is on shell number 4. The line segments pointing from shell n at time t to shell n + 1 at time 
t + At symbolize forward transfer, with magnitude proportional to the thickness of the line segment. Line segments pointing 
from shell n + 1 at time t to shell n at time t + At symbolizes backward transfer. 

Figure 2 shows the non-linear transfer, A„ = fc n _iim(un_iit n u,i+i) of energy from shell n— 1 to shells n and n + 1 
in the 3D case. The abscissa is the shell number and the ordinate is time, read from top to bottom, the figure is 
composed of line-segments connecting shells n at time U and shells n + 1 at time U + At, symbolizing an energy 
transfer from shell n in the time- interval [U ; ti + At] . The thickness of the line-segment is proportional to the size of 
the transfer. Line-segments going from n + 1 at time to n at time ti + At symbolizes the inverse energy transfer. 
The figure shows that the transfer is temporally intermittent and occurring in bursts B] , with the transfer being faster 
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as the bursts propagate to higher wave-number shells. From this the interpretation is straight forward; The residence 
time for a burst at a given shell is proportional to the eddy turnover time resulting in a more and more intermittent 
transfer as the eddy turnover time decreases. In order to illustrate how this kind of behavior can lead to anomalous 
scaling behavior, consider the following linear stochastic model, which is an extreme case. 

Let x n be a stochastic variable representing the energy at shell n, governed by the dynamical equation; 

x n (t+l) =p„_i(i)a;„_i(i) - p n (t)x n (t) - vk 2 n x n (t) + fS nA , (2) 

where p n is a stochastic variable which is 1 with probability tm /t„ and with probability 1 — tn jr n , thus it represents 
a transition probability of energy transfer from shell n to shell n + 1. The boundary conditions are X-\ = pn = 0. 
The average residence time is then simply proportional to the eddy turnover time. It is easily seen that x n is always 
positive. 
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FIG. 3. The result of a simulation of the simple stochastic model, x„ is plotted as a function of time. The intermittent 
transfer is similar to what is seen in the GOY model, figure 2. 

Figure 3 is similar to figure 2, but shows the values of x n . There is obviously no inverse transfer in this model. By 
comparing figures 2 and 3 one sees that the intermittent structure of the transfers are similar. The energy will in this 
model scale inversely proportional to the eddy turnover time, so that if we identify x n with the energy /enstrophy of 
the GOY model E n = £^K| 2 , the scaling of (\u n \) ~ fc^ (a+1)/3 is the same as for the GOY model. However, the 
structure function changes completely. It is readily calculated, with the result; (x^) ~ fci" 2 ^ 3 x q , where t„ — fci" 2 ^ 3 
is the eddy turnover time at shell n and x — fr\ is the average energy input into the system during one large eddy 
turnover time. This means that the structure function, ((q), becomes independent of q. q will only show up in the 
offset, glog(a;) in a log(x'), log(fc„) plot. This is illustrated in figure 4, showing the result of a numerical simulation. 
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FIG. 4. Anomalous scaling behavior of the simple stochastic model. log((:r 9 )) vs. log(fc) for q — 1,8. The spectral 
slope does not depend on q, So that £(9) = (a — 2)/3. The model was iterated for 10 4 large eddy turnover times with 
AT = 20,/ = 0.05,1/ = 10 -11 . 

Obviously this linear stochastic model cannot reproduce the scaling exponents of the GOY model or of real turbu- 
lence but it qualitatively illustrates the effect of temporal intermittency. 

In conclusion, the behavior of the 2d like GOY model shows either statistical equilibrium or cascade of enstrophy 
depending on the ratio of the eddy turnover timescales between the shells. The dynamics and the multiscaling of the 
3d like GOY model depends on the existence of the second non-positive definite inviscid invariant, helicity, through 
the way the helicity is dissipated in the model. In both cases, the multiscaling is an effect of temporal intermittency 
also originating from differences in eddy turnover times at the different scales, thus the 2d model with a = 2 shows 
no anomalous scaling. The effect is illustrated in a simple stochastic model. The scaling behavior of this simple 
model does not correspond to what is seen in the GOY model or in real turbulence. However, the model points to a 
mechanism of how the temporal intermittency can lead to anomalous scaling behavior. In order to further clarify the 
relationship between the spatial intermittency observed in turbulence and the temporal intermittency seen in these 
simple models, it would be interesting to see if the structures seen in the time - wave- number domain, figure 1, can 
be identified in direct numerical simulations of the Navier-Stokes equation. 
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